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Universal constants and 
equations of turbulent motion 

by Helmut Z. Baumert, 

Institute for Applied Marine and Limnic Studies (lAMARIS e.V.), Hamburg, Germany 

Abstract. In the spirit of Prandtl's [1926] conjecture, for turbulence at 
Re — 7- oo we present an analogy with the kinetic theory of gases, with dipoles 
made of quasi-rigid and "dressed" vortex tubes as frictionless, incompressible 
but deformable quasi-particles. Their movements are governed by Helmholtz' 
elementary vortex rules applied locally. A contact interaction or "collision" 
leads either to random scatter of a trajectory or to the formation of two 
likewise rotating, fundamentally unstable whirls forming a dissipative patch 
slowly rotating around its center of mass which is almost at rest. This 
approach predicts von Karman's constant as k = l/A/27r ~ 0.399 and 
the spatio-temporal dynamics of energy-containing time and length scales 
controlling turbulent mixing Baumert , 2009 in agreement with observations. 



A link to turbulence spectra was missing so far. In the present paper 
it is shown that the above image of random vortex-dipole movements is 
compatible with Kolmogorov's turbulence spectra if dissipative patches, 
beginning as two likewise rotating eddies, evolve locally into a space-filling 
bearing in the sense of Herrmann 1990 , i.e. into an "Apollonian gear" 



consisting of incompressible and flexibly deformable vortex tubes which are 
frictionless, excepting the dissipative scale of size zero. 

For steady and locally homogeneous conditions our approach predicts 
the pre-factor in the three-dimensional Eulerian wavenumber spectrum, 
E{k) = aie^/^ /c"^/^, as ai = |(4 7r)^/^ ~ 1.802, and in the Lagrangian 
frequency spectrum, E^u) = f3ieuj^'^, as (3i = 2. The unique values for 
«!, Pi and K are situated well within the broad scatter range of observational, 
experimental and approximative results. Our derivations rest on geometry, 
methods from many-particle physics, and on elementary conservation laws. 



1. Introduction 

In the present paper we show that fluid turbulence 
can be understood in an idealised sense as a statis- 
tical many-body ensemble - a tangle of vortex tubes 
taken as discrete particles. We follow an early conjec- 



ture by Ludwig Prandtl 1926 who discussed an anal- 
ogy between molecular diffusion and turbulence. He 



related his mixing lengtl[^ or Mischungsweg with the 
mean-free path of kinetic gas theory and considered his 
fluid elements or fluid lumps (Flussigkeitsballen) of lo- 
cally nearly same size as relatives of gas molecules. He 
further assumed his mixing length to scale with the "di- 



^According to Hinze 
the notion mixing length. 



1953 



also G. I. Taylor made early use of 
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ameter" of his fluid elements. 

Although purely heuristic, his concept became popu- 
lar in the years before WW2. The question how to com- 
pute the details could not be answered without reference 
to measurements. After WW2 the gas-kinetic analogy 
found thus strong criticism from the continuous-image 
side [see Introduction in Batchelor 1953 . This was 



plausible because the gas analogy represents a discrete 
concept that, if useful, would label a Copernicanian 
turn or a paradigm shift in our grown view of turbu- 
lence as an exclusive matter of continuum mechanics. 

in the foot- 



In this sitution Werner Albring 1981 



steps of Prandtl, explicitely challenged the continuous 



paradigm of Reynolds 1895 



1924 , Taylor 1935 , Batchelor 



Keller and Friedmann 



1953 and their many 



followers, when he posed the fundamental question: 

Can the Navier-Stokes equation be used to 
calculate turbulent flows? 

His doubts were based on the feeling that the continuous 
image of RANS does possibly override features of most 
elementary vortex interactions at small scales. 

In analogy to Albring one might ask: Can we deduce 
a rose's blossom from the periodic system of chemical 
elements? Intuitively we answer with no. However, this 
answer is justified because generally "more is different" 



Anderson 1972 



Looking in this sense at Albring's above question 
we may add that the prediction horizon of the Navier- 
Stokes equation (NSE) is strongly limited by a series 
of higher-order non-equilibrium phase transitions when 
a growing Reynolds number goes through a number of 
critical values. In these super-critical regions the sen- 
sitivity against initial conditions becomes relevant and 
leads to the famous butterfly effect Lorenz 1963 . This 



is closely related with irreversibility. At least in the 
limit of vanishing viscosity the initial- value problem for 
NSE looks like a reversible one. But we know that tur- 
bulence at Re = oo is an irreversible process. 

In principle this problem was already known by 
Poincare and others as an aspect of the three- 
body problem of celestial mechanics, but in a fluid- 
mechanical context it has been demonstrated only af- 
ter WW2 by ILorenzl 119601 119631. Without going into 



details of the onset of turbulence and non-equilibrium 
phase transitions in hydrodynamics and their reflections 
in various branches of pure mathematics, together with 
Prandtl and Albring we hypothesize that NSE is not 
sufficient to understand the secrets of turbulence or to 
even predict turbulent flows, in particular not in the 



limit Re — > oo. Some readers might even go a step 
further and argue that higher-order elements of the 
Friedman-Keller expansions of NSE, e.g. the second and 
third-order turbulence closures discussed in | Voropayeva| 



2007 and the literature quoted therein, as well as the 



various "corrections" of those closures populating tur- 
bulence theory, are modern analogues of the epicycles 
of geocentric times. But this would exceed the limits 
of the present report. In this respect we better refer 
to the agreeable picture of today's status of turbulence 
science drawn by Davidson 2004 in the preface of his 



book where he mentions even "religious wars. . . between 
the different camps" of turbulence theory. 

There are not many physical phenomena resisting 
theoreticians so long like turbulence. In the past it 
was mainly supra-conductivity which took about 30 
Today similarly obstinate prob- 



Feynman 1963 



years 

lems are dark matter, dark energy, and super-symmetry 
which have also now an age of about up to 40 years. But 
turbulence waits still much longer for redemption and 
remained particularly as the very last unsolved enigma 
of classical physics. 

Below we develop an asymptotically invariant alter- 
native to RANS and higher-order closures, which is not 
primarily based on NSE but does not violate NSE ei- 
ther. 

Although the theory of many-particle physics offers 
a huge reservoir of potentially helpful methods and 
tools, a closer inspection reveals that the number of 
directly applicable tools is less impressive. Neither Li- 
ouville theorem, ergodic hypothesis or Hamilton for- 
malism nor other concepts for thermodynamic equilib- 
rium are applicable in their classical forms. Turbu- 
lence is essentially an open-system, thcrmodymically 
non-equilibrium phenomenon. In best cases we have 
a steady state or Fliefigleichgewicht in the sense of 



Bertalanffy 1953 , Glansdorff and Prigogine 1971 and 



Haken 1978 1983 . However, the concepts of Ising 



about particle dressing, quasi-particles and renormal- 
ization - in the broader sense of Dresden 1993 and 



McComb 2004 



appeared as essential guidelines in 
the slow evolution of our thoughts. 

Surely, turbulence can successfully be attacked from 
more than one sider] Whereas the continuous image has 
led to a number of important results, for instance RANS 
and the broad spectrum of different heuristic closure 



^Theories generally compromise portrait and design aspects 
and in some case mathematically different images of turbulence 
will show up eventually as fully equivalent in their physical pre- 
dictions, like Heisenberg's matrix mechanics and Schrodinger's 
wave mechanics. 
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schemes discussed e.g. by Wilcox 2006 , without ad- 



ditional phenomenological input it gave neither unique 
values of the universal constants of turbulent motion 
like von-Karman's and Kolmogorov's spectral constants 
nor a closed image of turbulent flows. This has led to 



the following view shared by many theoreticians Lan- 
dauandLifshitz[|1987[ p. 173]: 



. . . and K (is) a numerical constant, the 
von Karman constant, whose value cannot 
be calculated theoretically and must be de- 
termined experimentally. It is found to be 
K = 0.4. 

In this respect one can actually be more optimistic. 
Frictionless turbulence is governed in a weak sense by 
the Euler equation, a special case of NSE. The Euler 
equation together with the conservation of mass repre- 
sent pure "inert geometry". We can thus expect that 
the universal constants of turbulent motion are ruled by 
certain geometrical constants, e.g. by irrational num- 
bers like v^ and/or transcendental numbers like tv. In- 
deed, this is possible. In a precursor study based on the 
particle picture used also below, von Karman' s constant 
could be derived as k = 1 /\/27r ~ 0.399 



Baumert 



2009 . This interesting result led to the conjecture that 



the particle concept can give also insight into spectral 
aspects of turbulence. 

Below we make explicit use of theoretico-physical 
thinking, which differs from other scientific activities 
like computer modelling and mathemtical physics, from 
applied mathematics, and from the measuring and ob- 
servational disciplines. 

Computer modelling and measurements have in com- 
mon that they deal with really existing, finite systems 
and data, in particular with real and integer numbers, 
with variables of finite size. 

Theoretical physics, however, deals with "things" 
which never did, which do not and never will exist: 
e.g. with point masses, homogeneous continua, ideal- 
ized vortex tubes, linear waves, frictionless fluids, plane 
and impenetrable walls, infinitely small or large vari- 
ables. A further building block of theoretical physics is 
an often unspoken principle which we apply also here: 
it "dictates that, all things being equal, one goes for the 
simplest possibility-a rule that has worked remarkably 
well." 



Zee 1986 



p. 85, see also Chandrasekhar, 1979]. 
Sometimes this principle is called Occam's razor. 

Mathematics also deals with similar non-existing ob- 
jects, but their relations to the real world are outside 
its responsibility. It is the responsibility of theoretical 
physics to relate idealized thought systems with the real 



world. The following report concentrates solely on the 
latter aspects. 

With one exception, the following text contains no 
"fancy shmancy mathematics" . It is the relation be- 
tween equations ([T]), Q and ([9| that would need some 
time to be derived from scratch. Furtunately, this has 
already been done by other authors many years ago and 
went into the textbooks on stochastic-dynamic systems 



e.g. Kraichnan 1968 Haken[ 1978 1983| Stratonovich 


1992| 1994 to which we refer the theoretically inter- 



ested reader. We further suggest to begin first with 
Baumert 2009 because it is actually the basis of the 



present report. 

Everywhere the pronomen "we" is used in this text, 
it means the two of us, the dear reader and the author. 

2. Particles 

If we take Prandtl's discrete particle concept seriously 
but not literally, then a number of questions need to be 
posed ~ and answered. A first question to be addressed 
is the following. 

Where do the particles come from? With re- 
spect to their substance our particles arc indistinguish- 
able from their fluid or gaseous environment. It is only 
their state of motion which makes the difference to the 
sourrounding fluid so that we better should talk about 
quasi-particles, but keep for brevity the term particles 
in the following. They are characterized by their ge- 
ometry, their specific kinetic energy, by their vectors of 
linear and angular momentum. 

Our particles are generated as representatives of tur- 
bulent fluctuations such that their energy has to be 
drawn from the mean flow field via so-called shear pro- 
ductiorjj to keep the total kinetic energy of the flow in 
balance. I.e. if there is a shear in the mean flow and 
non-vanishing turbulent viscosity, then the mean flow 
looses kinetic energy which reappears in form of TKErl 
i.e. as energy of vortex motions. Hereby an exact pro- 
cess description of vortex generation is not needed as 
long as we know how much energy is lost and how much 
vorticity is generated so that we can place vortices of 
corresponding properties into the flow. In a way we 
may talk about an emergence of the particles. 



^Convective turbulence and internal-wave breaking will be 
treated in another study. 

^In the following the word turbulent kinetic energy or TKE is 
often used. In the given context it means actually a kinetic energy 
density, i.e. energy per unit of mass. Therefore we measure TKE 
most comfortably in the units m^ s~^. 
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What do the particles look like? Due to 

Helmholtz' principle of conservation of circulation, a 
circulation-free fluid volume should exhibit zero circu- 
lation over the course of time. I.e. particle generation 
cannot add any circulation to the volume. This implies 
that turbulence production can take place only in form 
of the generation of couples of counter-rotating vortices 
with zero total circulation. Consequently we specify 
Prandtl's fluid elements as dipolcs where the simplest 
form is locally fully symmetric, which is to be under- 
stood in a statistical sense. 

This assumption is central for our theory but de- 
serves a comment. It is well-known that in laminar and 
turbulent flows wall friction generates velocity gradients 
or current shear, which is identical with circulaton in- 
duction into the mean flow. This process is governed by 
RANS, more speciflcally by its turbulent friction term 
containing correlators like {u' w'), i.e. by the Reynolds 
tensor. This tensor is calculated further below on the 
basis of a particle theory. If a particle would have own 
circulation on the local level, it would add extra cir- 
culation to the flow, which would violate RANS. Our 
particles can therefore have no own circulation. 

Thus we have the challenging theory situation that 
a mean flow exhibits turbulent friction and circulation 
while the turbulent vortices responsible for fluctuations 
and fricition have no circulation and are essentially fric- 
tionless, excepting the singular dissipation scale of size 
zero, see further below. 

With respect to the specific form of vortices in our 
dipoles we refer to the classical notion of vortex tubes 
wherein vorticity is confined to the interior of a tube 
of smaller or larger cross section [see further below and 



Pullin and Saffman 1998t who quote papers by Kuo 
& Corrsin, 1972, and Brown & Roshko, 1974, show- 
ing tubes as dominating characteristic structures]. A 
nice and more recent study of vortex tubes has been 
presented by Wilczek] [2011 who also published on the 
internet a number of valuable animations of vortex en- 
sembles in motiorj^ We note that, as always, the cen- 
terlines of our vortex tubes form either closed loops or 
they are attached to boundaries. 

Ordinary linear momentum of a dipole is known from 
classical theory in local approximation and can be imag- 
ined in analogy to the motions of smoke rings. It is 
conserved in our many-particle image if each particle, 
generated in a locally homogeneous volume element, ex- 
hibits no preferred direction of motion. The linear mo- 



mentums of two dipoles with opposite directions com- 
pensate each other. Assuming a locally "sufficiently 
high" and even number of dipoles thus guarantees that 
the particle-induced linear momentum and circulation 
are zero and total momentum of the flow exclusively 
governed by the mean-flow momentum balance includ- 
ing the Reynolds stress terms. 

How do the particles move? The motion of a 
single isolated dipole [a "naked" particle in the ter- 
minology of many-particle physics, see e.g. |Dresden[ 
McComb 2004 can be given by classical rules [see 
Lambl 119321 |A bring[ [T98T| |Saflman[ [T9921 [Baumert 



2005b, 2009 . However, if the selected dipole is embed- 
ded in a dense tangle and thus surrounded by a cloud 
of similar dipoles, its properties are "screened" by the 
cloud and thus modifled: it appears to be "dressed" - 
without violation of governing conservation laws. The 
solution is discussed in the next section. 



3. Vortices, energy, and scales 
Traditional vortex models 

The hydrodynamic literature presents a larger number 
of elementary analytic models for single isolated vor- 
tices under different conditions. Whereas the so-called 
potential vortex is less realistic, the Oseen vortex, the 
Rankinc and the Taylor vortex, the Burgers, the Lund- 
gren, and the Long, the Sullivan and the spherical Hill 
vortex are more realistic models of isolated vortices 
far from others and from boundaries [for overviews see 



LaSIbl |T932l [Higtl [19791 |Albring[ \mT\ |Saffman[ [1992 



Puflin and Saffman 1998 Davidson 2004 



^http : //pauli . uni-muenster . de/tp/menu/f orschen/ 
ag-f riedrich/mitarbeiter/wilczek-michael.html 



Excepting the potential vortex, all other vortex mod- 
els are principally realistic and show first a core with ra- 
dially increasing tangential velocity, then a saddle, and 
then a tail wherein the tangential velocity decreases ra- 
dially down to zero. However, in our context these mod- 
els are not applicable because they hold for conditions of 
isolation only. They cannot be transferred to conditions 
of a dense vortex tangle where the distance between vor- 
tex dipoles is small and the surrounding cloud of similar 
vortices screens the effects from the rest. Therefore a 
new approach is needed, without violating governing 
conservation laws, neither on a global nor on the scale 
of the locally homogeneous and isotropic fluid volume. 

I.e. an ideally "dressed" vortex dipole moving fric- 
tionless within a tangle of similar objects should nec- 
essarily be characterized by a finite effective radius, r, 
within which all the kinetic energy, /C, and vorticity 
of the vortex is concentrated so that we may talk of an 
"energy-containing radius" . For simplicity, in a statisti- 
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cal sense r and uj are taken identical for the two vortices 
forming the internally symmetric dipole. 
We summarize as follows: 

• The effective scales uj and r are governed by local 
conservation of energy and angular momentum. 

• The radius r defines a boundary within which ki- 
netic energy and vorticity of a vortex are confined. 

• The effective tangential velocity of a vortex is 
u = Lur and equals the propagation velocity of the 
dipole. The resulting kinetic energy of the dipole 

is thus /C = 2 X u'^/2 = 2 x r'^uj'^/2 = (rw)^. 

• The vortices are incompressible but deformable 
quasi-solid bodies which move frictionless in the 
vortex tangle. I.e. vorticity is uniformly dis- 
tributed within the cross-section of a vortex tube 
and has the value 2 x w. 

• At solid boundaries vortices perform frictionless 
roll motions. Dissipation happens exclusively 
within the fluid at scale zero. (In reality the dis- 
sipation interacts with boundaries through heat 
and sound generation. 

These assumptions describe our image of vortices in a 
turbulent vortex tangle without violating the conser- 
vation laws of kinetic energy, momentum, angular mo- 
mentum and circulation. Compared with the above- 
mentioned vortex models our approach may be inter- 
preted as a renormalization procedure leading to a finite 
kinetic energy and a finite spatial extension. 

While on a first glance this procedure looks some- 
what arbitrary, exactly this image is well established 
since long in physical oceanography, which is a world of 
highest Reynolds numbers. In physical oceanography, 
meteorology and physical limnology the a.m. character- 
istic radius, r, is traditionally related with the "energy- 
containing" scale, which is a direct observable under 
conditions of stable stratification and related with the 
so-called overturning and Thorpe scales. Central as- 
pects are discussed in detail in the following section. 

Vortex tubes and rigid-body rotation 

In stable stratification turbulent vortices are directly 
recognizable in vertical profiles of density, p, and other 
scalar variables, such as temperature and salinity, as 
"overturning" of density profiles. Within certain en- 
ergetic limits they move heavier over lighter fiuid and 
thus produce a direct imprint of the vortices in instan- 
taneous depth profiles of density, p{z). This imprint 



has a unique character in stratified fiows because the 
buoyancy force implies that there should be a stably 
stratified "reference" state, p, a state in the absence of 
turbulent overturning. 

This imprint appears typically as a Z-pattern, see 
Fig. [T] which obviously stands for a solid-body rota- 
tion wherein the tangential (azimutal) velocity grows 
linearly with the radius coordinate up to a maximum 
at its outer radius and then (outside the body) drops 
rapidly down to zero. We interprete a Z-pattern as a 
realisation of a vortex embedded in a dense vortex en- 
semb l e. i.e . as a so-called rigid vortex as discussed by 
Lugtl (19791 and further belowfl 



To our knowledge, the use of overturning deviations 
from the stable state to derive vortex radii was first 



introduced and explored by Thorpe 1977 in an anal 



ysis of temperature profiles from Loch Ness, Scotland. 
Thorpe showed how the reference state can be recon- 
structed from measurements and how overturns and 
overturning scales can be detected. His approach is a 
solid basis to measure the radii (or diameters) of vor- 
tices in stably stratified fluids directly, 
described 



Thorpe 1977 



the general case of a 
monotonous reference density distribution and the over- 
turning fluid body as frozen for the time of the over- 
turn. This is justified as long as the length achieved 
through molecular diffusion during the time of rota- 
tion. Id = \/2 Vrn tot 1 is small compared with the energy- 
containing radius r of the overturrTl In water this con- 
dition is always satisfied as Id remains here in the range 
of centimeters and less. 



Later authors like Imberger and Boashash 1986| con- 



sidered a subset of the monotonous case, the linear 
reference-density distribution as depicted in our Fig. [l] 
Itsweire et al. 1986 in their Fig. 1 were the first to 



discribe the overturning motion explicitely as the effect 



of a rigid-vortex motion [cf. also Pullin and Saffman 
19981. 



Thorpe's method in oceanography 

Background. A short digression on the background of 
Thorpe's method might be in place. Geophysical fluid 
dynamics is a domain of pioneering turbulence stud- 



® The idea of a rigid vortex is compatible with fluid mechanics 
because Helmholtz' laws also apply here. But the connection with 
continuum mechanics is not trivial because the velocity distribu- 
tion in a rigid vortex exhibits a singularity at the outer radius 
where it drops from a finite value down to zero. Further, the 
isolated rigid vortex is not stable. 

^Here Urn is molecular viscosity and tot the duration of an 
overturn at the energy-containing scale r oc L. 
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Figure 1. Overturning and overturning scales in the 



sense of Imberger and Boasliasti 1986, Fig. 6 tlicre 



The cartoon depicts a 180° overturning of a 20 m thick 
layer that brings relatively heavy water up and lighter 
water down. Note the characteristic Z-pattern of the 
vertical displacement (. The original stable density pro- 
file is the blue line, the unstable density in the overturn 
is red and dashed. By courtesy of Dr. Hartmut Peters, 
Earth and Space Research, Seattle, USA. 



ies at high Reynolds numbers because atmosphere and 
oceans offer the necessary conditions for free. The first 
massive demonstration of measured Kolmogorov spec- 



tra was carried out by Grant et al. 1959 at a Reynolds 
number of about Re « 10** in a 100 m deep tidal chan- 
neP] They demonstrated the validity of the 5/3 law 
over an intervall of 3.5 orders of magnitude (about 12 
octaves). Together with the logarithmic law of the wall 
and the decay of TKE in homogeneous isotropic turbu- 
lence according to t~^ [or {x/U)~^ in the wind tunnel; 
of. Batchelor 1953 Chapter VII] the 5/3 law belongs 
until now to the most prominent universal features of 
high-i?e turbulence. 

Many observational techniques and dynamical con- 
cepts of geophysical fluid dynamics were first developed 
in meteorology and later adopted in physical oceanog- 
raphy and limnology. That this was not the case with 
turbulent overturning is likely related to the fact that 
the troposphere, the most accessible part of the atmo- 



^Today's DNS of turbulent flows do not exceed Re Ri 10^ and 
even the European superpipe CiCLOPE will in th e best case not 
significantly exceed Re ft: 10^ Riiedi et al. 2009 . 



sphere, does not show consistent mean stratification. 
In contrast, oceans and lakes are almost everywhere 
stably stratified outside of boundary layers. It fur- 
ther makes sense that turbulent overturning was first 
explored in a lake where temperature, T, is the sole 
stratifying agent. In contrast, ocean stratification de- 
pends on temperature as well as salinity, S, and den- 
sity p = p{p, T, S) cannot be measured with the same 
resolution as T alone. We note in passing that real flu- 
ids are compressible such that analyzes of overturning 
scales have to be based on potential density and poten- 
tial temperature. 

The imprint of turbulence in vertical density profiles 
allows defining and extracting overturning scales such 
as the Thorpe scale, Lth Thorpe 1977 . The cartoon 



of Fig. [T] illustrates the overturning of a Lot = 20 m 
thick layer of the water column by a single vortex in 
solid body rotation of diameter Lot ■ According to our 
above discussion of vortex dipoles as quasi-particles we 
conclude that Lot = 2 r. 

The graph depicts the moment after a 180° rotation 
that brings heavy water up by vertical displacements 
with a range of C = ... -I- Lot and moves lighter water 
down by C = ... — Lot- The Thorpe scale is defined as 
the r.m.s. of C. The linear original density profile of the 
cartoon implies 



L 



Th 



LotlVi 



Vi' 



1.15 r 



(1) 



We define C, by the path from its original depth to a 
displaced depth. That is, C, carries the same sign as the 
vertical turbulent velocity w' that has caused it. 

Fig. IT] may seem simplistic. It can easily be made 
more realistic by adding that, in the course of the over- 
turning, and owing to the unstable stratification within 
overturns, flow instability will occur and generate tur- 
bulence and a range of smaller scales than Lot, inside 
the big overturn. 

The ocean is full of overturns that look just like 
this scenario. Fig. [2] depicts a big overturn in the Pa- 
cific Equatorial Undercurrent, EUC, on the equator at 
140° W from IPeters et all 119951. Note the Z-shape of 



the big overturn and its sharp upper and lower edges. 

In oceanography and limnology, Thorpe's concept of 
vortices overturning parts of the water column are ap- 
plied to observations as in Fig. [2] Measured poten- 
tial density (ae; by convention 1000 kg m^"^ is sub- 
tracted) or potential temperature (O) data points are 
"Thorpe-sorted" into monotonically rising or falling se- 
quences corresponding to stable density stratification. 
The sorted profiles are taken as a proxy for the refer- 
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Figure 2. A big overturn in the Pacific Equatorial 
Undercurrent at 0°,140°W during the Tropic Heat II 
cruise [Peters et al.[ 1995 adapted from ]: (a) poten- 
tial temperature (O, red) and potential density (tre, 
black), Thorpe-sorted ce (blue), and Thorpe scale Lrh 
(shaded); (b) turbulent vertical displacement (. By 
courtesy of Dr. Hartmut Peters. 



ence, or "mean" profile that gave rise to the observed 
turbulence and overturning. The vertical distance over 
which data points have to be moved to make the pro- 
file monotonic is —(. I.e. Thorpe-sorting "undoes" 
the overturning. Individual overturns are defined by 



Dillon 1982 



The shaded bars in Fig.[2ja) show Lxh averaged over 
individual overturns. This definition is highly robust. 
A squared buoyancy frequency iV^ computed from the 
sorted ae or Q is non-negativq^ N"^ > 0. 

Turbulent length scales related to overturning scales 
are unaffected by the presence of internal gravity waves 
(IGWs). IGWs are ubiquitous in stratified geophysical 
flows and dominate velocity and scalar spectra at ver- 
tical scales of the order of 0(100 m) and smaller [e.g. 
[Peters et al. 1995 . The TKE accounts only for a tiny 
fraction of the integral of these "red" spectra. Owing 
to this property, length scales commonly used in lab- 



^ Here we may restrict our considerations to incompressible 
fluids where N'^ = —g/{p)d{p)/dz with g as gravitational accel- 
eration, g = 9, 81 m s~^. 



oratory experiments, such as the Ellison 1957 scale, 
Le = p' /{d{p) /dz) , are unsuitable for a characteriza- 
tion of turbulent length scales in geophysical flows be- 
cause density fluctuations p' = p — (p) inevitably are 
dominated by IGW signals so that Le can no longer be 
interpreted as a characteristic length scale of turbulent 
fluctuations. Here, angular brackets denote ensemble 
averages p' — {p'^)^^^ is an r.m.s. value. 

Thorpe scale and Taylor scaling. The great 
power of Thorpe's [1977] concept can be demonstrated 
by considering the so-called Taylor scaling, with q as 
an r.m.s. measure of turbulent velocity fluctuations. 
Dimensional analysis strictly gives for I as an energy- 
related length scale the following: 



I ex q^/e 



(2) 



Here I is obviously related with TKE because ofJC oc q"^, 
and £ is the TKE dissipation ratq^ 

If the Thorpe scale is indeed a measure of the size 
of the energy-containing vortices and thus related with 
the energy-containing scale, i.e. Lth oc /, then we can 
write ([2]) as follows. 



£ X Lxh oc q"^ 



(3) 



g ex (£ X LthV^^ . (4) 

Instead of q one can alternatively use K}'^ with tur- 
bulent kinetic or mean vortex kinetic energy, K,. Fig. 
[3| shows measurements in the ocean demonstrating 
that relation Q is true, in a statistical sense. The 
Thorpe scale characterizes indeed the size of the energy- 
containing eddies. 

However, this statement needs a comment. Thorpe- 
sorting as a technique to measure the radius or diame- 
ter of energy-containing turbulent vortices rests on the 
existence, and hence works only, in stratified flows. Un- 
der geophysical conditions this automatically implies 
that a certain coexistence of IGWs and turbulence is 
inescapable. The g in (pi) must only reflect TKE and 
must not be contaminated by the much larger IGW en- 
ergy. 

While the separation of waves and turbulence is 
beyond the scope of this notq^ a workaround is to 
consider spectra and to study only the ultra-violett 



'^"We use the convention that the symbol /C denotes the energy 
of one selected dipole whereas K. is the local average of the same 
variable for a whole ensemble 

'^^The reader is referred to 
[Lien|[2000[ . 



Peters et al. 



1995 



and 



D'Asaro and 
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Figure 3. Turbulent velocity qf — {eiy^^ derived from 
measured e via Taylor scaling, (HI), versus measured tur- 
bulent velocity fluctuation q. Adapted from |Peters et al.| 
1995 where l = l.QLth and where q is the spectral ve- 
locity variance at vertical wavenumbers > l/l. Shown 
are data for individual overturns at least 1 m thick 
from depths 60-350 m in the Equatorial Undercurrent 
at 0°, 140° W. The gray line indicates the median ratio 
of ge over q. By courtesy of Dr. Hartmut Peters. 



or short-wavelengths part where IGW existence is ex- 



cluded through their dispersion relation. [Peters et al 
1995 did exactly this and extracted q from oceanic 



observations as the velocity variance at scales of I and 
smaller. On this basis, Taylor scaling clearly emerges 
from their observations (Fig. Isl) . 

Thorpe's [1977] concept of turbulent vortices and 
their imprint through vertical overturning on density 
profiles allows extracting energy-scale variables even 
in conditions where velocity and scalar spectra are 
heavily dominated by internal gravity waves. It al- 
lows relating energy-scale- to dissipation scale variables 
through the turbulent length scale that carries Thorpe's 
name. These directly observable energy-containing 
length scales stand for the radius or diameter of our 
renormalized vortices discussed further above. 

Dressing a rigid-vortex tube 

As we learned from the considerations and examples 
above, turbulence may be imagined as vortex tubes re- 
sembling a dense "local cloud" of dipoles in chaotic mo- 



tion, the cloud having zero angular and linear momen- 
tum. So far the rigid-vortex tubes have been consid- 
ered as somewhat arbitrary idealizations of real-world 
vortices. 

Translated into the language of home cooking, a 
snapshot of rigid- vortex turbulence may be imagined as 
a dense, entangled heap of hot spaghetti arrabiata. The 
spaghetti rotate around their inner centerlines and move 
frictionless within the (inviscid) sauce arrabiata. There- 
fore the only interactions between individual spaghetti 
occurs when they touch each other randomly. 

It arises the question whether the rigid-vortex tube 
and a corresponding dipole is a stable solution of the 
Euler equations. We leave it to the interested reader to 
derive from these equations that an isolated rigid vortex 
indeed solves these equations, but, as long as it exists, 
generates the following pressure head as a consequence 
of inertial (centrifugal) forces: 



P^Po 



P 2 2 

- X uj'^r'^. 
2 



(5) 



Here po is the background pressure of a laminar refer- 
ence flow, e.g. in the ocean the depth-depending hydro- 
static pressure. If the pressure outside the vortex would 
simply be pg then, due to the action of the outwards- 
directed pressure head of the vortex motion given by 
the second term in ([s]) , the vortex would loose stability. 

The stability of quasi-rigid vortices as observed in 
real-world turbulence can thus be guaranteed only by 
the help of a compensating force of equal strength. Here 
a concept of many-particle physics comes into play: 
dressing. As far as we will embed our initially isolated 
rigid vortex into a locally homogeneous and isotropic 
large ensemble ( "cloud" ) of similar vortices (more pre- 
cisely: dipols made of rigid vortices), the members of 
the cloud generate more or less exactly the counter pres- 
sure needed to compensate ([5]) and keep the vortex "sta- 
ble enough" . 

According to ([5| and with u — cor the pressure de- 
viation can be written as follows. 



Sp = p-po 



P 2 

2 



p X JC, 



(6) 



where /C is the kinetic energy density in a vortex. Ob- 
viously the measurement of turbulent pressure fluctua- 
tions can help estimating turbulent kinetic energy. This 
corresponds to an early classical but approximate result 
of continuum theory [eq. (8.3.21) on p. 182 in Batch- 

242, last eq.]. We will not 



elor 1953 Hinze 1953 



explore the potential of (l6| further and leave it to the 



interested reader. 



Baumert: Universal constants and equations 



4. Turbulence and kinetic gas theory 
Similarities 

About 150 years ago, James Clerk Maxwell presented 
his kinetic theory of gases. Even if the details of molecu- 
lar interaction forces in vacuum would have been known 
that time, it would have not been helpful. An inte- 
gration of Newtons law of motion in terms of ordinary 
differential equations for each one of the billion parti- 
cles was practically excluded. But the creative use of 
symmetries and conservation laws radically simplified 
the situation and led eventually to a closed description 
of the most important macroscopic properties of gases. 
Main elements are the following: 

(i) The particles in a gas are perfectly elastic points 
with non-zero mass. 

(ii) They arc in permanent random motion which is to 
be described in terms of statistical moments and 
sometimes called "molecular chaos". The scatter 
motions exhibit no preferred dirctions. 

(iii) The collision results depend only on the local an- 
gular orientation of the collision partners. 

(iv) Between collisions they move uniformly and inde- 
pendently, without preferred direction. 

(v) Due to chemical neutrality, collisions lead only to 
the scatter of trajectories. 

The corresponding "turbulent relatives" of the above 
are: 

(i') The particles are locally symmetric vortex-dipole 
tubes with finite cross-sectional area, with vortic- 
ity and kinetic energy confined in the tubes, and 
with zero circulation. 

(ii') Their random translatory motions prefer no direc- 



tions and may be termed 'dipol chaos' Marmanis 



1998 



(iii') The result of collisions depends only on the lo- 
cal angular orientation of the colliding dipole ele- 
ments. 

(iv') Between collisions the particles move along com- 
plex trajectories which may be curved. 

(v') For symmetry reasons, 50 % of all collisions oc- 
cur between two counter-rotating dipole elements 
leading to dipole recombinations [or reconnec- 
tions, like those reported in a turbulent super- 



quasi-elastic, random scatter motions resembling 
turbulent diffusion and mixing, see left branch in 
Figure |4] 

We underline that these rules shall apply only locally in 
space and time. Of course, the prediction of the global 
pathway of a vortex dipole is impossible. 

Differences 

Major differences between ideal gases and turbulence 
are the following: 

(a) The number of chemically neutral gas particles in 
a countainer is constant in time. They form a 
closed thermodynamic system. There is no parti- 
cle annihilation. 

(b) The particles of turbulence are excited energy 
states, their number decreases via collisions and 
subsequent energy dissipation. They need to be 
replaced by shear generation of new particles if 
their number shall be kept more or less in bal- 
ance. 

(c) In a gas, trajectories of point masses are simple 
lines in space, trajectories of vortex-dipole tubes 
produce curved areas. 

While gases evolve towards static thermodynamic equi- 
libria, turbulence evolves either towards dynamic equi- 
libria or steady state^^ which may have periodic char- 
acter, or turbulence dies off. 

To specify this important difference we have to sup- 
plement our above property list for "turbulence parti- 
cles" as follows: 

(d) For the same symmetry reasons like in (v') of sec- 
tion [4] above, the remaining 50 % of all dipole col- 
lisions occur between two likewise rotating dipole 
elements. Each case generates a fundamentally 
unstable vortex couple forming a slowly rotating 
dissipative patch with its center of mass more or 
less at rest. This patch decays through turbulent 
dissipation in the special form of a "devil's gear" 
or Kolmogorov spectrum. For details see further 
below and the right branch in Fig. [3] 



^^This neither implies homogeneity, equifinality nor uniqueness 
of steady states. 



fluid by Paoletti et al. 2010 and at the end to 
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free dipoles with 
zero total circulation: 







interaction 



*N 



\"»» "*-/ ^-k- »» »T 



recombined stable dipoles: 
turbulent diffusion 



Instable configurations: 
turbulent dissipation 



Figure 4. Two possible collision results of two dipoles: 
the left branch is "diffusive" because it leads to a recom- 
bination of the dipole elements and scatter of trajecto- 
ries which is known as turbulent diffusion. The right 
branch is "dissipative" because it leads to an unsta- 
ble vortex configuration which decays "somehow" into 
heatp^ For symmetry reasons both branches have iden- 
tical probabilities of 0.5. Note that the circular form of 
the vortex cross sections presented here are chosen for 
reasons of clarity. They hold as an statistical average 
picture only. Real vortices have elliptic or even strongly 
deformatted cross sections. 

5. Equations of turbulent motion 
Brownian and turbulent motions 

The turbulence properties (iii') and (v') of section H 
above establish an analogy of turbulent dipole move- 
ments with the Brownian motion of particles suspended 
in a fluid at rest [in the sense of Einstein 1905 . In the 



present turbulent case it is not the kinetic theory of 
heat which governs the motions, it rather is Helmholtz' 
theory of dipol motions. 

The (local) temporal path increments SYj of a vortex 
dipole j may be found by integration of the following 
Langevin equation over the time increment 6t, 



dt ' ' 



J = 1...AA, 



(7) 



where Vj is the random center-of-mass velocity of the 
selected dipole j. 



To specify the stochastic process Vj = Vj (t) as simple 
as possible, we choose a zero-mean, white-noise Gaus- 
sian process. Its strength is controlled by the locally 
averaged dipole properties /C and Co. More elaborate 
random processes like the Ornstein-Uhlcnbeck or the 
Kraichnan model Kraichnan 1968 and others are be- 



yond the focus of the present paper. 

We now need to make a bigger jump over the broad 
river of stochastic dynamic systems theory where the 
probability density function for the solution of a dy- 
namic systerrrjmay easily be taken from applied text- 
books like those of Haken 1978 1983| in form of solu- 
tions of the Fokker-Planck equation (FPE) correspond- 
ing to the extremely simple stochastic-dynamic system 
([t]). In the present case the FPE reads as follows. 



dl\f d 
dt dx 



dx 



= 0, 



(8) 



where Af is the probability density mentioned above, an 
equivalent of the particle number density itself, i/ char- 
acterizes the strength or intensity of the noise Vj (i) in 
([t]) . As already mentioned it is governed by the control 
variables (/C, ui) and appears later as coefficient of eddy 
diffusivity of momentum, i.e. as eddy viscosity. 

In the transition from (It]) to (Is]) it has implicitely 
provided that all gradients exhibit sufBciently smooth 
and slow behavior. Those local quasi-equilibrium con- 
ditions are typically assumed in non-equilibrium ther- 
modynamics and many-particle physics. This clearly 
excludes shocks and steep fronts from our considera- 
tions. 

A further comment concerns an assumption used im- 
plicitely above. It is an analogue of the so-called Ising 
assumption known from the theory of magnetism: we 
treat triple or higher interactions between vortex fila- 
ments as negligiblq^ 



^^A dynamic system is here understood as a set of possibly 
non-linear ordinary differential equations driven at their right- 
hand sides by stochastic processes. 

^^This is not fully trivial. Under specific conditions of Bose- 
Einstein condensates stable configurations consisting of one vor- 
tex and two anti-vortices have been observed in the laboratory, 
either in linear setups or equilateral triangles [Seman et al.[[2009| . 
The latter is most symmetric and called a tripole. Here we as- 
sume that in our very dense "vortex gas" tripole-tripole interac- 
tions are controlled by dipole-dipole interactions of their subsets. 
Note that tripoles might violate the principle of zero circulation 
on the level of elementary interactions needed to keep the mean- 
flow momentum balance in correct balance. 
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Figure 5. Local cross section through the first de- 
velopmental stage of an unstable pair of likewise ro- 
tating vortices resulting from a dipole-dipole collision 
(right branch in Fig. [4]). The green circles represent 
the primary energy-containing vortices with identical 
radii r. Note that the green circles do not touch each 
other! They are separated by the red circles who sym- 
bolize secondary vortices in the beginning phase of a 
whole vortex cascade. The dipole cloud surrounding 
the above structure not only generates the necessary 
pressure head to keep individual vortices stable but it 
also acts as a source of perturbations initiating roll-up 
instabilities and thus tertiary and higher-order vortices 
and eventually a fully developed dissipative patch, see 
text. The broken blue line and the arrows symbolize 
the slow rotation of the whole patch around the com- 
mon center of mass. 

Generation and annihilation of particles 

Property (d) in Section H above allows to supplement 
the right-hand side of equation (Is]) with sink and source 
terms, leading to a reaction-diffusion type equation: 



dAf 
~dt 



d_ 
dx 



m 

dx 



= n - /3 A/'2 



(9) 



Here 11 is the rate of quasi-particle generation. It is to 
be expressed in terms of kinetic energy per unit time 
and is thus necessarily proportional to the energy loss 
of the mean flow. 



The second term corresponds to the energy- 
dissipation rate of TKE, e, i.e. TKE conversion into 
heat and/or sound, with /? being a constant [for details 
e scales with the rate of particle 



see 



Baumert 2009 



annihilation resulting from collisions. As we learn from 
chemical kinetics [see Haken 1978 1983 for details], 



the annihilation term is quadratic in the particle num- 
ber, Af, because two particles need to collide to generate 
one unstable couple which is then converted into heat 
and/or sound. This is discussed in the next Section. 

6. Dissipative patches 
Formation of vortex spectra 

Our Fig. [4] and the concept behind it may seem sim- 
plistic. It demonstrates the only two possible results 
of a dipole collision. While the left half of the Figure 
shows the recombination of cownter-rotating vortices 
from counter-rotating vortices, the right half shows the 
the formation of a couple of likewise rotating vortices 
from counter-rotating vortices. If isolated or naked, the 
likewise rotating couple revolves around a common cen- 
ter of mass and remais thus nearly at rest. This couple 
is known to be unstable, to form stationary dissipative 
19481^ ^1 



patches Sommerfeld 



This picture with the dissipative patch almost at 
rest implies that dissipation is a spatially patchy phe- 
nomenon called intermittency and studied extensively 



by various authors [for an overview see e.g. Frisch 1995 



We do not go into details because intermittency is out- 
side the focus of this study. We discuss instead the way 
how the unstable configurations at the right half of Fig. 
Hand in Fig.lslcould be transformed into a Kolmogorov- 
Richardson spectrum. This problem has attracted early 
attention by |Taylor| |l937| and |Kolmogorov| [l94l] . The 
latter found on dimensional grounds that for a steady 
energy flux from large to small scales the kinetic energy 
spectrum as function of wavenumber may be presented 
as follows^ 



dK. = ai 6°'^ k- 



dk 



(10) 



where here k — 2ir /\ is the wave number and A the 
wavelength, e is the dissipation rate of TKE, K.. Based 
on strict dimensional arguments, [Kolmogorov 1941 
proved that a2 = 2/3 and as = 5/3, in excellent 



^^It somehow resembles the "dissipative elements" discussed by 
the group around Norbert Peters [seelSchafer et al.||2010|. 

'^^Actually, the original arguments in Kolmogorov's derivation 
were more subtle, but the use of arguments in the sense of 
Rayleigh's method of dimensional analysis or, stricter, Bucking- 
ham's TT theorem, is "strict enought" for the present discussion. 
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agreement with the famous observations by [Grant et aL| 



1959 in a tidal inlet with Re ~ 10 and a depth of 
about 100 m. A theoretically sound value of ai was 
open until today and given now below. 

Devil's gear 

Our view of the Kolmogorov-Richardson cascade has 



been changed through a study by Herrmann 1990 who 



has shown that Kolmogorov's value for a^ corresponds 
numerically to the data of a space-filling bearing [see 
also Herrmann et al. 1990] . The latter is the dens- 



est non-overlapping (Apollonian) circle packing in the 
plane, with side condition that the circles are pointwise 
in contact but able to rotate freely, without friction or 
slipping. One may call it a "devil's gear" Poppe 2004 



The contact condition for two different "wheels" with 
indizes 1 and 2 of the gear reads 



UJi ri = LJ2 r2 



(11) 



where u is necessarily constant throughout the gear and 
governed by the energy of the decaying vortex pair as 
u = y/2IC. It follows that 



UJ2 



ri 

UJi — 

r2 



(12) 



so that for very small r2 the frequency uj2 may become 
acustically relevant. 

If the above gear is frictionless then the next ques- 
tion arises where - within this picture - energy could 
be dissipated. Clearly, in a fluid with vanishing but 
non-zero viscosity, dissipation happens at scales where 
velocity gradients are high enough, here: at a scale of 
measure zero. Our dissipative patch (Fig. [s] shows the 
first stage of its formation) is thus "almost frictionless" 
and therefore a Hamiltonian system, excepting scales of 
size zero. 

The formation of a fully developed spectrum of 
"wheels" from Fig. [5] deserves certain perturbations 
"from the sides" , a condition which is guaranteed by the 
random reconncction/recombination and scatter pro- 
cesses sketched in the left half of Fig. [2] and also by the 
incomplete mutual pressure compensation of the vor- 
tices in our vortex ensemble. 

Without speculating too much we may expect that 
in a quasi-steady state patches like in Fig. [s] are formed 
via roll-up instabilities at the "surface" of the respec- 
tive larger vortices. They steadily evolve into a fully 
developed gear. Its energy, dissipated at the smallest 
radii, will decrease the energy content of the primary 
(initial) vortex pair unless it is fed by mean flow. 



The outer limits of such a patch are sketched in Fig. 
[6]for the begin of the cascade process. The most impor- 
tant message which we gain from this figure is that the 
longest or energy-containing wavelength of the dissipa- 
tive patch equals Aq = 2 r. The wavelength in a dipole 
is 4 r and it forms no patch or spectrum. This difference 
is essential. Further below we use Aq as a lower inte- 
gration limit for the spectral energy distribution. It is 
important to underline that Aq labels the upper wave- 
length limit in a dissipative patch. This limit is actually 
not influenced by the formation details of the spectrum. 

We finally notice that the transformation of the un- 
stable configuration of two likewise rotating vortices 
deserves time to set the greater masses of the smaller 
scales into motion. This inertia effect might play a role 
in highly dynamic scenarios. 




Figure 6. Outer limits of a dissipative patch (c./. Fig. 
[5]). The maximum wavelength is obviously equal to 
Ao-2r. 



7. Universal constants 



In a precursor study Baumert 2009 it has been shown 



that the geometric-mechanical concept given in the 
present paper allows the derivation of von-Karman's 
constant as k = I/VStt. It fueled hopes that other 
universal constants of turbulent motion might also be 
derived from that apparatus. 

Kolmogorov constant ai 

in the wavenumber spectrum 



The TKE can be calculated by integrating ( 10 ) over the 



dissipatve patch in the sense sketched in Figures [5] and 
[6] yielding 



/C 



ai e 



2/3 



fc-5/3 ^fc 



ai 



2/3 



(13) 



where fcp = 2 tt/Aq characterizes the lower end of the 
turbulence spectrum in the wavenumber space. We 
loosely assign the wavenumber range k ~ . . . ko to the 
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mean flow which may basically be resolved in numerical 
models. 



The dissipation rate e in ( 13 1 can be expressed as 
follows, 

e^lC/r, (14) 

with T being the lifetime of a dissipative patch. Insert- 



ing (14) in (13) and rearranging gives the following 

2 ,^ ,9/^..w, / r \2/3 
3 



ai 



(2.f/^/CV3(_L)^ 



(15) 



Here we took from Fig. [6] that the energy-containing 
initial or primary wavelength of a dissipative patch is 
given by Aq = 2 r. 

Now we make a local quasi-equilibrium assumption 
for conditions of extremely dense vortex packing: Our 
marching dipoles can occupy only those places which 
are simultaneously "emptied" from dissipative patches 
by decay. This means that the life time of a dissipative 
patch, r = /C/e, should equal the time of "free flight" 
of a dipole over a distance 2r [cf. Baumert 2009| : 

r = /C/£ = 2r/u. (16) 



Here we used the scalar dipole velocity u, 



: ujr 



(17) 



After some algebra, the dimensionless pre-factor of the 
three-dimensional wavenumber spectrum reads as fol- 
lows: 

Q-i = ^(47r)2/3 = 1.802. (18) 

The corresponding value of an ideal one-dimensional 
spectrum is one third of the above, i.e. 0.60. 

Kolmogorov constant /3i 
in the frequency spectrum 

The generalized form of the Lagrangian frequency spec- 
trum of fluid turbulence is the following. 



dK. = /3i e^ 



-/33 



duj , 



where uj is the angular frequency. |Tennekes and Lumley| 
derived 6 



1972 



1 and /33 = 2 in a similar fashion 



like Kolmogorov 1941 derived a2 and as [see also Me 



Comb 2004 . The integration of (19) from uj = ojq to 



UJ 



oo gives after rearrangement and using ( 14 ) the 



following. 



/?! 



Wo T . 



(20) 



Now we take from the right part of ( 16 ) t — 2 r /u, from 



the left part of (17) wg — ujr and insert both in (20) 
to get finally 

u 2r 



Wo r 



2. 



(21) 



Pre-factor in the velocity 
autocorrelation function 

The spatial autocorrelation function B{p) of fiuctuating 
velocities is a special second-order case of a structure 
function^ With the abbreviation p — |p|, it is defined 
as follows. 



B{p) = (Ui(x) X M2(f +/?)). 



(22) 



Here u\ , u^ are the (scalar) velocity components of the 
flow velocity u in the direction p connecting the points 
X and X -\r p where the velocities ui and U2 are taken 
respectively: 



Ml 
U2 



u{x) ■ p/p , 
u{x + p) ■ p/p . 



(23) 
(24) 



The central dot denotes the scalar product (or dot prod- 
uct). Notice that, rather than a density, p is here a 
spatial distance. 



In their §34 on p. 145 Landau and Lifshitz 1987 have 



shown that, based on early results by Kolmogorov, B{r) 
in ( 22 ) may be written as follows. 



B{r) ^ C X (s pf/^ 



(25) 



where C is a dimensionless numerical constant which is 
related with ai from the universal wavenumber spec- 



trum ( 10 ) as follows. 



c = «i|r(i/3) 
= ^(4^)^/3 1 r(i/3) 

w 2.37 . 



Here T{z) is the Euler gamma function. 

In contrast to our derivations of k, ai and /3i 



(26) 

(27) 
(28) 

the 



(19) value (28) for C should be taken with care because the 



derivation of its relation with ai by Landau and Lifshitz 
uses an approximation and is thus valid for small values 
of the distance variable p only, i.e. for ^ Ao = 2 r. 

8. Discussion 

Comparison virith observations 

The rounded numerical values k = 0.4, ai = 1.8, 
ai/3 = 0.6, and /?i = 2 for von Karman's and Kol- 
mogorov's universal constants predicted by our theory 
are situated well within the error bars of many high-i?e 



For details see 534 in Landau and Lifshitz 1987 
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number observations, NSE and RG based analytical ap- 
proximations, laboratory and DNS experiments. Based 



of the estimate of /3i seem unlikely. 



on observations, Tennekes and Lumley 1972 gave the 
values ai = 1.62 and /3i = 2.02, but with greater un- 
certainty. 

Later works have been analysed in an important 



study by Sreenivasan 1995 who possibly gave the most 



comprehensive literature review of experimental and ob- 
servational values for the number ai/3 until now. Later 
Yeung and Zhoul 119971 reported a value of ai = 1.62 



based on high-resolution DNS studies with up to 512^ 
grid points. 




Figure 7. Experimental and observational results 
for ai/3 measured, collected from the literature, and 
analysed by iSreenivasanl 11995 . The solid green line 



follows our somewhat arbitrary approximation 0.6 x 
\/Re\/ (^s/Rel -\- \/Re\) wherein ai/3 = 0.6 is the the- 
oretically derived asymptotic value. Here we took 
-Re* = 10. A similar presentation has been chosen by 



Sreenivasan 1995 for ai/3 in his Fig. 3, and by Yeung 



et al. 2006 for 7r/3i in their Fig. 2. 



The results of their later efforts suggest, again with 
DNS but based on a grid of 2048^ points, the value 
/?! = 2.1 Yeung et al. 2006 their Fig. 3]. In a recent 
study by Donzis and Sreenivasan 2010] a DNS grid of 
4096^ has led to ai « 1.58. 

Based on their oceanic measurements (with much 



higher Reynolds numbers compared with DNS) Lien 
[and D'Asarol [2002] found that ^i = 1.75 . . . 2.04. They 
state that 

. . . since the present uncertainty is com- 
parable to that between high quality esti- 
mates of the Eulerian one-dimensional lon- 
gitudinal Kolmogorov constant measured by 
many dozen investigators over the last 50 
years, large improvements in the accuracy 



Beginning with an initiating work by [Forster et al.| 
[1977 , systematic analytical approximations using RG 



methods and related techniques for NSE became fur- 
ther sources of estimates for the universal constants. 
E.g. Yakhot and Orszag 1986b ^ reported ai w 1.62 
whereas McComb and Watt 1992 derived ai = 1.60 ± 
0.01 and Park and Deem 2003 obtained ai = 1.68. 



We note that these approximations are technically ex- 
tremely complex and neither unique nor part of an in- 
tegrated descriptive concept for turbulence. 

Turbulence plays a crucial role in almost all fields of 
engineering, including medical applications, and in geo- 
physical fluid dynamics up to climate-change studies. 
It plays an essential part in our everyday life. When 
we leave our house or ride our bike, when we jump into 
our pool - we always are literally embedded in a turbu- 
lent fluid. Compared with most other fields of modern 
physics the present scatter in the values of the univer- 
sal constants is uniquely high and therefore actually not 
longer acceptable. However, our theoretical results con- 
taining irrational and transcendental numbers suggest 
that the universal constants of turbulence can princi- 
pally not be measured in real-world fluids. Or would 
measuring tt be a reasonable task? 

A possible solution of this dilemma is the asymp- 
totic analysis of measurements at higher and higher 
Reynolds numbers so that an extrapolation to Re — > oo 
becomes feasible with some certainty. The data analy- 



2006 



ses by|Sreenivasan 1995 and the report by Yeung et al. 
could serve as a methodical model. Fig. [7] shows 



a re-plot of Sreenivasan's data with Re on the abscissa 
and measured ai/3 on the ordinate. The data exhibit a 
visible tendency to grow with increasing Re to a satura- 
tion value which is statistically indistinguishable from 
our 0.6. 

Universality and fundamentality; 
turbulence, physics, and geometry 

Besides problems and tasks posed in engineering and 
geophysical fluid dynamics where prediction accuracy 
and extrapolability matter most, there is a question 
posed by natural philosophy: how are universal con- 
stants of turbulent motion, fundamental physics con- 
stants like the electron's elementary charge [see e.g. 
and mathematical constants like tt re- 



Fritzsch 2009 



lated with each other? 

First, the measurement accuracy of "turbulent num- 
bers" is extremely poor. 

Second, "turbulent constants" characterize universal 
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properties of a specific (turbulent) form of dynamic and 
self-similar motions rather than more "static" proper- 
ties like the elementary charge. They are thus closer to 
the vacuum speed of light or the Hubble constant de- 
scribing the (accelerated) expansion motion of the uni- 
verse. 

Third, like mathematical constants, "turbulent num- 
bers" are dimensionless whereas fundamental physics 
motion constants like the Hubble or the vacuum speed 
of light are given in kilometers per second. The lat- 
ter are absolute values. "Our" constants characterize 
self-similar motional 

This provocates the question whether geometry is 
part of physics or vice versa, as discussed for instance 
by Palais 1981 . One might argue that the discovery of 



where (p) is the pressure. Equation (31 ) is usually called 



the constant angle sum in triangles was historically the 
first discovery of a physical conservation law by man, 
made and explicitly formulated much earlier than the 
conservation laws of volume, mass etcp^ 

The completed image of turbulent flov^rs 

For the sake of clarity we concentrate on the most sim- 
ple non-trivial situation, a spatially one-dimensional 
channel flow with velocity component U in horizontal 
(x) direction, with vertical variation along z, and with 
horizontal, u, and vertical velocity fluctuations, w. An 
example of this situation has been presented earlier [see 



eq. (4.20) and (4.21) in Baumert 2005a where stratifi 



cation is already considered through the squared buoy- 
ancy frequency, iV^]. The Reynolds decomposition of 
our flow field then reads as follows: 



U{z,t) = {U)+il{z,t), 
W{z,t) = {W)+w{z,t). 



(29) 
(30) 



Mean flow: RANS. One may now insert ([29) [30]) 
into the corresponding two-dimensional Euler equation 
to find together with mass conservation the following. 



d{U) d{iiw) 
dt dz 



d{p) 
dx 



(31) 



^^ If we exclude from our consideration physical cosmology and 
arbitrarily chosen ratios between masses of the various atoms and 
molecules then physics has actually only one dimensionless excep- 
tion: Sommerfeld's fine-structure constant, with a value of about 
137. Wolfgang Pauli has long been preoccupied with the ques- 
tion of why, and Richard |Feynman| [1985| even speculated about 
a relation with tt. Here we have shown that tt is at least related 
with constants of turbulent motion. 

■^^ These thoughts are actually outside the scope of the present 
report. Some patience is needed as a "theory of everything" is not 
in sight [Laughlin] |2005| , and just-answered questions typically 
give birth to new conundrums. 



a Reynolds- averaged Navier-Stokes equation (RANS). 

Here we take the case of vanishing molecular viscosity. 

Turbulent mixing: downgradient flux. Equa- 



tion (31) is not yet closed because the correlator de- 
scribing a diffusive flux, {Hw), needs to be specifled, a 
task which consumed substantial efforts over the last 
60 years. The following quasi-linear flux-gradient rela- 
tion is theoretically well established in many branches of 
many-particle physics and reads in our case as follows: 



[uw) 



oz 



(32) 



Kolmogorov-Prandtl relation. However, also 



( 32 1 does not yet close the problem because now the 
so-called turbulent viscosity, v, needs to be specifled. 
In our picture the latter is again presented in a many- 
particle format in analogy to Einstein's theory of Brow- 



nian motion as follows [for details see Baumert 2009|: 



L^/T = JC/{7rn) 



(33) 



Here L — t/k and t = l/fl are locally averageq | space 
and time scales which are expressed within the frame- 
work of our mechanistic dipol-chaos model in terms of 
/C and w = 27r f2. 

Equations of turbulent motion, neutral strat- 
iflcation. We skip here the technical details of deriva- 

and quote only the result 
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tions given in 

that the variables /C and fi are governed by a specific 
system of nonlinear partial differential equations as fol- 
lows: 



die 


d 


dt 


dz 


an 


d 


dt 


dz 



dK. 



dfl 



m 



n' 



2\ dz ) 



(34) 



(35) 



These equations resemble at least the structure of the 
so-called K,-Vt closure model used by many authors [loc. 
They use differing empirically gained 



Wilcox 2006 



cit 

sets of prefactors of the terms and cannot identify the 



physical nature of Vt. Our derivation of (34 351 should 



not be confused with a new "scheme" or closure of this 
sort because - as we have shown - our equations are 
founded solidly on the most simple principles of hy- 
drodynamics and many-particle physics, without use of 
phenomenological data and, for the first time, giving 
even the universal constants of turbulent motion. On 



For reasons of transparency the overbars are omitted. 
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the other hand, the experience of Wilcox 2006 and col 



leagues shows that practical experiments and applica- 
tion needs guide the creative engineer very close to the 
physically correct solution, of course, without guiding 
further to the universal constants. 

Turbulence spectra. With the dissipation rate e, 

e^lCn/TT, (36) 

we may express the spectra in the admissible turbulent 
ranges k = kg . . . oo and lo — luq .. .oo as follows: 

die 



dk - "1^ ^ 


(37) 


d-^ o -2 

-— = (3ieuj ^ , 
duj 


(38) 


turbulence begins at 




fco = tt/kL = TT^n^/2/IC , 


(39) 


UJQ = n/K^ = 27rJ7, 


(40) 



respectively. Longer scales or slower motions are possi- 
bly turbulent, but not in the sense of a fully developed 
spectrum. In the case of slow variations of the driv- 
ing forces the motions outside the validity range belong 
to the mean flow. Notice that the lower limits of the 
turbulent spectral ranges are dynamic quantities which 
are controlled by the dynamic variables K. and fl. 

Universal constants. The values of our theoret- 
ically derived universal constants of turbulent motion 
are 

1 

3 

2, 

/2^ = 0.399 . 



Pi 



(47r)2/3 = 1.802, 



1/^ 



(41) 
(42) 
(43) 



Solid boundaries. TKE cannot penetrate solid 
boundaries like walls, which is sometimes called an 'adi- 
abatic boundary condition', in analogy to heat conduc- 
tion. Assuming the wall at z = z* this means that the 
following condition has to be satisfied: 



die 
dz 



= 0. 



(44) 



This is already sufficient to solve (34 35 ) at solid walls 



and produce the logarithmic law of the wall. A second 
condition for f2 is not needed due to the a.m. nonlinear- 



ities of the equations [for details see Baumert 2005a 



Stable stratification. Stable stratification may be 



described by the following modification of ( 34 1 



dK d 
dt dz 



dK. 

'Ih 



d{U) \ 
dz ) 



2N^ - n^ 



(45) 



This description is valid as long a.s il > N. In the case 
of free decay and linear stratification, fl approaches N 
"from above" . Due to the dispersion relation for inter- 
nal waves, slow disordered turbulent motions at fl ^ N 
and slower can only exist as waves so that turbulence is 
converted into internal waves when J7 — >■ TV. The strat- 
ification aspects are discussed by [Baumert and Peters] 



2004 2005 in greater detail and in relation to obser- 



vations and measurements. 



Limitations. A trivial limitation for (45 1 is N > 0. 



For unstable or convective situations these equations 
are possibly not applicable, but we did not yet test this 
case. 



Equations ( 29 - 44 ) form a closed, complete descrip- 



tion of mean and turbulent motions, including their 
wavenumber and frequency spectra with their universal 
spectral constants, applicable under neutrally stratified 



conditions. For stably stratified flows (34) needs to be 



replaced with ( 45 ) , provided that there are no external 



sources of internal- wave energy. This means that, due 
to the ubiquituous presence of internal waves in most 



geophysical flows, the latter are not covered by (45 1 



Some modifications are necessary and a certain knowl- 
edge of the sources of internal-wave energy is needed. 
The necessary modifications are not part of this report. 



The difference between (34 35 1 or ( 45 35 ) on the one 



hand, and traditional phenomenologically based closure 



schemes discussed e.g. by Wilcox 2006 on the other, is 



mathematically minimal but physically relevant. Those 
schemes use for instance quantities like e or r as pri- 
mary variables in balance equations although conser- 
vation laws for those quantities do not exist. The dif- 
ferences are becoming generally relevant and visible at 
solid boundaries and, in the stratified case, also within 
the turbulent fiuid volume. 

An important limitation has already been mentioned 
above. Further, variations in the mean flow flcld should 
be relatively slow and spatial gradients not too strong 
so that we may talk about time-dependent but quasi- 
steady behavior with sufficiently homogeneous state 
variables on local scales and sufficiently developed Kol- 
mogorov spectra. This excludes shock-waves from our 
considerations. But alos this limitation has not been 
tested yet. 

The free decay of turbulence might be such a case. 
Our theory predicts the decay of TKE to scale with 
t~^ [or with {x/U)~^ in the wind tunnel]. This is 
also supported through a group-theoretical analysis by 



Oberlack 2002 . Whereas wind-tunnel data Batchelor 



1953 and selected free-decay measurements Dickey and 



Mellor 1980 agree surprisingly well with our theory, in 
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some experiments if rapid free decay the exponent is 
somewhat greater than unity. This situation is chal- 
lenging and needs a further study. 
We may summarize as follows: 

• The present theory is based on conservation laws 
and geometry only, without use of empirical data. 

• Besides spectrally integrated parameters like tur- 
bulent viscosity, dissipation rate etc., it gives full 
turbulent spectra, even under non-stationary con- 
ditions. 

• The theory predicts fundamental constants of tur- 
bulent motion. 

• The approach is limited to slow and smooth mean 
flows under neutral or stable stratification. 

• Internal gravity waves are taken into account only 
so far as they are intrinsically coupled with local 
current shear. 

Our methodology differs from past concepts. The latter 
use specific series expansions of the Navier-Stokes equa- 
tion into a system of partial differential equations for 
higher and higher moments in the perturbations, with 
the zero*'' hierarchy element being the Reynolds equa- 
tion (RANS). The k-e, the Mellor-Yamada, fc-r, the k-uj 
and the many other turbulence-closure schemes popu- 



lating the literature are examples [e.g. Baumert et al. 



2005 Wilcox 2006 



The present picture includes a perturbation step, 
too, but it stops at the zero*'' level, at RANS. Here 
the turbulent fluxes still remain unknown. We then 
took into account that each fluid element experiences 
an effective total force field consisting additively of the 
external mean field Ti acting in the total fluid volume, 
and of the field H controlled by nearest neighboring 
vortex dipoles ("Ising assumption"). 



n = n + H. 



(46) 



Our results show that H is locally best described by a 
simple white-noise stochastic force. This approach has 
successfully been used about 100 years ago by Langevin, 
Smoluchowsky and Einstein in theoretical analyses of 
Brownian motion. This assumption directly delivers 
Fokker-Planck equations for the expectation-value dis- 
tribution functions of particle number and of particle 
properties /C and fi. 

The theory applies exclusively to locally homoge- 
neous, locally isotropic and weakly unsteady turbulent 



flows. To extreme non-stationarities and/or sharp spa- 
tial gradients like in Shockwaves it is possibly not ap- 
plicable. As a rule, temporal changes of the mean flow 
should take place on time scales sufficiently greater than 



l/rj because otherwise the Kolmogorov spectra (37 38) 
are not yet well enough established. 

Coda. The theory completed here will hopefully not 
form a Procrustean bed in Saffman's (1977) sense [loc. 
p. 107] who even feared that "in 



Davidson 2004 



cit 

searching for a theory of turbulence, perhaps we are 

looking for a chimera ...". This view has just recently 



been enforced by Hunt 2011 as follows: "But there are 



good reasons why the answer to the big question that 
Landau and Batchelor raised about whether there is a 
general theory of turbulence is probably 'no'. " 

As the above smacks of defeatism, we pass the word 
to Sir Winston Churchill who demanded in a critical 
situation "We shall never surrender!" He surely would 
encourage a new generation of men and women to fol- 
low their own fresh ideas towards a general theory of 
turbulence at asymptotically high Reynolds numbers. 
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